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ABSTRACT 

Two different algorithms — a deterministic magnetic-field-only algorithm and a 
Kalman filter for gyroless spacecraft — are used to estimate the attitude and rates 
of the Rossi X-Ray Timing Explorer (RXTE) using only measurements from a 
three-axis magnetometer. The performance of these algorithms is examined 
using in-flight data from various scenarios. In particular, significant 
enhancements in accuracies are observed when the telemetered magnetometer 
data are accurately calibrated using a recently developed calibration algorithm. 
Interesting features observed in these studies of the inertial-pointing RXTE 
include a remarkable sensitivity of the filter to the numerical values of the noise 
parameters and relatively long convergence time spans. By analogy, the accuracy 
of the deterministic scheme is noticeably lower as a result of reduced rates of 
change of the body-fixed geomagnetic field. Preliminary results show the filter- 
per-axis attitude accuracies ranging between 0.1 and 0.5 deg and rate accuracies 
between 0.001 deg/sec and 0.005 deg./sec, whereas the deterministic method 
needs a more sophisticated techniques for smoothing time derivatives of the 
measured geomagnetic field to clearly distinguish both attitude and rate solutions 
from the numerical noise. Also included is a new theoretical development in the 
deterministic algorithm: the transformation of a transcendental equation in the 
original theory into an S^-order polynomial equation. It is shown that this 8 - 
order polynomial reduces to quadratic equations in the two limiting cases— 
infinitely high wheel momentum, and constant rates-discussed in previous 
publications. 


INTRODUCTION 

It has been demonstrated 1 " 10 that the attitude and rates of low-Earth orbiting spacecraft can be 

simultaneously estimated using measurements of the Earth's magnetic field, B , using only a three-axis 
magnetometer (TAM) and no a priori information. The feasibility of this "TAM-Only" scheme essentially 
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relies upon 2? changing direction rapidly enough in the spacecraft body frame to make computation of its 
time derivative possible, and these changes during the course of an orbit are sufficiently large to enable 
determination of all three Euler angles using only TAM data. 

Our approach consists of using two independent algorithms-deterministic attitude determination 
from magnetometer-only data (DADMOD) and the Real Time Sequential Filter (RTSF). The DADMOD 1 * 3 
is a TAM-only algorithm that relates the time derivatives of B in inertial and spacecraft body coordinates 
to determine the attitude and the body rates. The RTSF 4 ’ 5 is a robust Kalman filter that estimates, in 
addition to the attitude, errors in rates propagated via Euler's equation. Note that the RTSF is a general 
algorithm for gyroless spacecraft; however, its sensitivity to rate errors as small as 0.0003 deg/sec makes it 
a robust and accurate real-time algorithm even in TAM-only situations with no a priori spacecraft 
information. 

The highlights of our past applications to in-flight data from the Solar, Anomalous, and 
Magnetospheric Explorer (SAMPEX) and the Earth Radiation Budget Satellite (ERBS) have shown that 
with a T AM-Only approach: (1) SAMPEX attitude and rate requirements can be met even when the on- 
board Sun sensor fails 4-6 , (2) using partially calibrated magnetometer data from ERBS nominal mission 
mode the RTSF yields 8 accuracies within 0.4 deg and rate accuracies within 0.005 deg/sec, and (3) ERBS 
attitude and rates could be reliably determined 10 during its 1987 control anomaly 11 when the spacecraft 
tumbled at approximately 2 deg/sec. Another useful aspect of the past work is the combination of the 
strengths of these algorithms in an automated scheme' wherein the deterministic algorithm is used to 
initialize the more accurate Kalman filter to within a few degrees of the correct spacecraft attitude. 

In the present work, we examine the performances of these algorithms during some important 
scenarios of the RXTE: (1) calibrated and uncalibrated TAMs and (2) during maneuvers. While the first 
scenario does not require further explanation, the motivation for the second is the possible application of the 
RTSF to extend aging missions. For example, ERBS (launch: 1984) needs monthly thruster-based 
maneuvers for solar-power purposes, and these are currently conducted using rate information from the one 
remaining gyro channel. This paper demonstrates that, by providing magnetometer-only rate solutions, the 
filter can be a useful tool during such maneuvers, especially when the last gyro channel also fails. 

The present work concentrates on results for RXTE that are interesting in their own right, because, 
in contrast to SAMPEX and ERBS: 

• RXTE is inertial pointing so that B changes very slowly in the body frame, and this leads to 

observability and convergence issues when only short data spans are used (as is the case here). 

• RXTE is a zero-momentum spacecraft, whereas SAMPEX and ERBS are momentum-biased about 

the pitch axis. We believe this leads to the RTSF results being very sensitive to the numerical 

values of its propagation noise parameters and their relationship to the weightage of the TAM data. 

For a similar reason, the accuracy of the deterministic scheme becomes noticeably lower. 

The rest of the paper is organized as follows. Section 2 briefly describes the algorithms; included 
here are recent developments in DADMOD and a novel TAM calibration algorithm 12 that is used to 
calibrate the RXTE TAM data. Sections 3 and 4 analyze the performances for RXTE and ERBS 
respectively, and Section 5 summarizes the conclusions. 

ALGORITHMS 

Deterministic Attitude and Rate Determination Using Magnetometer-Only Data (DADMOD) 

As discussed in detail in previous publications, determination of the spacecra; attitude and rates 
based on magnetometer measurements and their first and second time derivatives can be cast in the form of 
the following vector equation: 
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( 1 ) 


A 0 (<J>) + A ,(0)©!+ A 2 a)! 2 = 0 

-A 

where the angle <t> and the body rate o>i around the body-fixed geomagnetic field vector# are two 
unknown variables. It is essential that 

= £ A x//,(<t>), A 2 =Q 2 x B a , (2) 

where the vectors H j (O) and Q 2 are given by Eqs. (5) and (6) in Ref. 10, and B A = B A / I B A I . 

The third vector A 0 (<I>) is also perpendicular to B A at any value of the angle <!> so that two nontrivial 

equations to determine both <I> and are obtained by projecting vector equation (1) on two directions 
perpendicular to the geomagnetic field. As a new development, we show that these two directions can be 
chosen in such a way that £(<D) = tan(<I>/2) becomes one of roots of an S^-order polynomial P 8 [Q . 

*A T 

In fact, by projecting Eq. (1) onto the vector B x A 2 one finds 

©,(<!>) = - I fi A l (A 0 (<D).O 2 )/(Hj(d>).i 2 ) (3) 


Substituting the latter expression into the projection of Eq. (1) onto the vector H j (O) then gives 

(H l (<t>) • A 2 ) ( A 0 (0>) . (d>)) + I B a I 2 (A 0 (O) .o 2 ) 2 = 0 (4) 

It can be shown that the vectors A 0 (<l>) and H l («I>) have the form: 

A 0 (<D) = c 4 (d >)0 4 [C(<P)] , H x (<J>) = c 2 (d>)G 2 [C(<I>)] (5) 


where the vectors Q a \Q are formed by polynomials of the n* order in £ and c(<I>) = cos (<l>/2) , and 
therefore the solution sought for is given by one of roots of the S^-order polynomial: 


where 

p 2 [q = 2 )* P4ta = i i (GJQ* & 2 ). Pfita s (G 2 K] *04^1) 


Note that all coefficients of polynomial (6) are equal to zero when the geomagnetic field is 
directed along one of the spacecraft principal axes of inertia since both vectors £2 2 and A 2 vanish. By 

analogy these coefficients are nullified if the vector A 2 becomes perpendicular to the vector H x (< t>,) for 
the sought -for root £(4>i). However, it can be shown that the ratio (P^O ) 2 / P 2 IX] tends to zero in both 
cases so that the solution sought for can be found among real roots of the polynomial P$[Q. Because Eq. 
(3) is no longer applicable, one has to solve a quadratic equation to find ©1 . To avoid instabilities, a 
special algorithm was developed to select the direction associated with the maximum of the discriminant. 
As a result, the resultant solution remains stable as the coefficient of the quadratic term tends to zero. 
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The ratio (P 4 [£] ) 2 / Pzt^] also vanishes as the wheel momentum h tends to while the 
polynomial P 6 [Q takes the form: 

p 6 [Q= d + ; 2 ) 2 n 2 tg (7) 

As a result, we come to the quadratic equation ri 2 [£] = 0 discussed in Ref. 6. A similar 
decomposition of the polynomial P 6 [£] takes place in case of constant body rates after one drops all terms 
associated with the time derivative of the angular velocity vector CO . The solution solved for can be found 

from the requirement for the vectors H j (<J>) and A 0 (<I>) to be perpendicular to each other, which is 

equivalent to the condition = 0. The resultant quadratic equation = 0 has been studied in 

detail in Refs. 1 and 2. 


Real Time Sequential Filter (RTSF) 

In view of space considerations, only details relevant to the tuning of the RTSF are presented here. 
A full mathematical description of the RTSF has been provided elsewhere (References 4 and 5). 

The RTSF’s state vector X is comprised of the four components of the attitude quaternion, q , and the 

corrections, b , to the spacecraft’s rates, a i : 

x = [r s r ] r < 8 ) 

(Note that the components of b and © are resolved along the spacecraft’s x, y, z axes.) 


The RTSF uses sensor data to estimate q as well as b , with b being estimated kinematically in 
the same manner as gyro biases for a gyro-based spacecraft; i.e., by attributing differences between the 


measured and propagated attitudes to errors in © . The b estimates are then used to correct © , and these 
corrected rates are used as initial conditions to propagate Euler’s equation to the next measurement time. 
The propagation of b is modeled via a first-order Markov model: 


db_ 

dt 


b - 


(9) 


where i) b is a white noise term, and t is a finite time constant. A suitable value for z is the time between 
measurements. 


The rates are assumed to contain a white noise component, fj a and are propagated using Euler’s 
equation after accounting for the angular momentum contributed by the wheels, and for the total external 
torques acting on the spacecraft. TAMONLY currently models the gravity-gradient torque and the magnetic 
control torque acting on the spacecraft. (The aerodynamic drag torque and the radiation pressure torque 
have been intentionally omitted to reduce the amount of spacecraft modeling required. The RTSF relies on 
the rate-corrections, b , to compensate for the small effects of these torques.) 

The covariance matrix, P, is propagated by numerically integrating the following equation: 

— = F(a>)P + PF T (a>)+Q (10) 

dt 
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Here F(co) is described in Reference 5; the quantity of interest is the 6 x 6 matrix Q that 
quantifies the propagation noise and is of the following diagonal form: 

Q = diag [ Q a , Q a , Q a , Q b , Q b , Qb ] ( l *) 

Here Q a is related to the noise term fj a and contributes to the growth of the attitude error 
covariances about the body X-, Y-, and Z-axes during propagation. Similarly, Q b is related to noise term 
Tj b , and contributes to the growth of the error covariances of b during propagation. Another quantity that 
we must consider during tuning is a, the strength of the white noise in the TAM measurements. 

The filter can be initialized in one of the following two ways before processing a span of telemetry data. 

• Inertial initial conditions (IIC), where the spacecraft is assumed at rest in the Geocentric Inertial 
Coordinates (GCI) with its axes coinciding with the GCI axes; this results in large initial errors. 

• Deterministic initial conditions (DIC) where the filter makes short (2 to 5 min) runs and determines 
which of the DADMOD solutions is a good a priori solution. This results in small initial errors. 

The TAM Calibration Algorithm 

The effects of TAM calibration were determined using a recently-developed algorithm 12 where the 
following set of 21 time-independent parameters are used to "adjust*' the magnetic field vector measured by 
the TAM whose axes nominally coincide with the spacecraft body axes. 


E 3x3 scale factor/misalignment matrix nominally equal to the identity matrix, / 3x3 

G 3x3 TAM-torquer coupling matrix nominally equal to the null matrix, 0 3x 3 

/ 3x1 bias vector nominally equal to the null vector, 0 3x 1 


If at any instant!? is the magnetic field vector measured by the TAM, D is the 3x1 dipole moment vector 
of the magnetic torquer bars, A is the known GCI -to-spacecraft body frame attitude matrix, and Bf is the 
corresponding Intemationl Geomagnetic Reference Field (IGRF) 13 vector in the inertial frame, the 
calibration model assumes 

B, = 1 (E 9 B*-G'Dj)-f t +v, , i= 1, 2, 3 (12) 

) = 1 

where B R = ABf is the predicted field in the spacecraft body frame, and o is a white-noise term of root- 
mean-square (r-m-s) value a. The goal of the calibration then is to estimate E, G, and / , by applying 
statistical methods to a span of TAM measurements, { B lt ...,B N } and the corresponding predictions 
{B?,...,B* }. 


Resolving E and G into the vectors E l , E 2 , E 2 , G, , G 2 , and G 3 as follows. 


E = 


w 

El 


G - 


Gj 


(13) 


three independent loss functions are now formulated as: 
L, =-^1 a Hj [Ej b; - Gj D.-fi- B ni f 

(T T ; fl=t L 


where the subscript n denotes the measurement time. 


i= 1, 2, 3 


(14) 
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(15) 


1 " 1 

1 Z 2 

.-i cr. 


and 


<2„ , =* 


# r.» 


The following notation for the statistical quantities formed from vectors is followed here. 
Means: = £ CL i X i 

' ' i»l 

Covariances: ^ f'j = (XY T ) ~ (x^(f T ^ 

where the superscript T denotes matrix transpose. 

Minimizing L^Ej, G, , / j yields 


(16) 

(17) 

(18) 

(19a) 

(19b) 
(19c) 

where subscript i indicates that the averages are only over the i-th set. Equations (19a) and (19b) can be 
readily solved yielding: 

or -[(Wrr -( s .| 6 ),(®"l®)''I(¥) 1 (*'l">r • -(o|«'), (*'| b' < 20 ») 

fi is then obtained by using Equations (20) in (19c). 

RESULTS USING THE ROSSI X-RAY TIMING EXPLORER IN-FLIGHT DATA 
Overview of the Mission and Data 




The RXTE is an inertial-pointing spacecraft and was launched in December 1995 into a near- 
circular orbit of altitude 580 km and inclination 23 deg. The primary attitude sensors on board are charge- 
coupled device star trackers that provide accurate sensor-determined attitudes during inertial periods. The 
attitude during maneuvers (as many as eight each day) is obtained from accurately calibrated gyros. The 
predicted field values B were generated using a 10 th order IGRF model for the reference field values. 
Three sets of data from 1/4/96, 7/4/96, and 1 1/6/97 were used in the present study. Of these, the first two 
contain spacecraft slews (primarily about the z axis), while the last is wholly inertial. The telemetry data 
received at the FDF are nominally 2 sec apart, but various samples at a slower rate were generated to 
increase observability of the magnetic field variations. Thus, data were generated with pseudo-periods 
ranging from 4 sec to 40 sec, and several different telemetry periods were used for each set of data. 
However, the results presented here used 40 sec sampling for the 7/4/96 data and 8 sec sampling for the 
other two sets. 


Terminology 


Some notes about the figures and tables presented here are in order. The "truth" models used to 
evaluate the attitude and rate accuracies of the algorithms are the on-board computer (OBC) determined 
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attitudes and rates computed from their time derivatives. GCI-to-Body attitude results are presented in the 
form of 1-2-3 Euler angles, and these angles are respectively referred to as "Angle- 1", " Angle-2", and " 
Angle-3". The body-frame components of the spacecraft rates are depicted in the figures as “wx", "wy", and 
"wz". "Raw" and "adjusted" refer to the quality of the TAM data, and denote pre- and post-calibration 
values for the TAM measurements. "Residuals" are the differences between TAM measurements and fields 
predicted using the RTSF attitude estimates. "TAM angle" is the angle between the measured and predicted 
fields. It is a convenient scalar parametrization of the separate TAM residuals along the three body axes 
and, as will be seen, is useful when evaluating the filter in the absence of truth models. Only TAM-1 
measurements have been used throughout the paper although TAM-2 measurements are also available for 
the RXTE. The TAM-2 measurements and the residual statistics are not very different from the TAM-1 
measurements, although significant differences do exist in the calibration parameters. "RTSF rate-errors ” 
are the corrections, b , estimated by the RTSF (see Equation (9)) as part of its state vector and are different 
from a term such as "error in wz" that refers to the differences between the RTSF rates and gyro rates. Thus 
b may be viewed as "rate residuals" since convergence of the RTSF implies small b . 


TAM-1 Calibration Results 


Excellent residual statistics were obtained after calibration of the data and the results are shown in 
Table 1 for each axis separately. For example, the root-mean-square residuals are of the order of 0.5 mG. 
The mean residuals are most impressive: of the order of 10' 14 (i. e. of the order of "e-14" in the notation of 
the Table). 


Table 1 

RXTE RESIDUAL STATISTICS FOR TAM-1 


Dataset 

Pre-Calibration (Raw) Data 
(mG) 

Post-Calibration (Adjusted) Data 
(mG) 


Mean Residuals 
(X,Y,Z) 

R-M-S Residuals 
(X,Y,Z) 

Mean Residuals 
(X,Y,Z) 

■ 

1/4/96 

4.105,-0.481,-1.792 
Max: 9.290 
Min: -6.949 

4.720, 1.563, 
3.191 

-2e-14, -9e-15, -6e-14 
Max: 3.019 

Min: -1.296 

0.519, 0.335, 
0,272 

7/4/96 

3.055,-4.510, 3.413 
Max: 7.984 
Min: -3.196 

3.162, 4.54, 
4.114 

2e-14, 2e-14, 4e-14 
Max: 0.712 

Min: -0.923 

0.139, 0.216, 
0.203 

11/06/97 

-1.561, -0.579,8.221 
Max: 15.670 
Min: -6.635 

2.778, 2.160, 
9.085 

7e-14,7e-14, 2e-13 
Max: 2.799 
Min: -1.544 

0.601,0.312, 

0.494 


RTSF Tuning 

The RTSF was tuned as follows. The largest of the r-m-s residuals results for a given dataset of 
Table 1 was used as the RTSF tuning parameter a during the TAM-only runs; for example, this value would 
be 4.72 mG for the raw data of 1/4/96. At the outset of the TAM-only runs approximate numerical values 
for the filter propagation noise parameters, Q a and Qb of Equation (11), were obtained by analyzing the 
errors in the angular momentum of the spacecraft and wheels and the effects these errors would have on the 
RTSF rate and attitude while propagating between measurements. The uncertainties in the system net wheel 
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angular momentum was determined to be about 0.025 N-m, which resulted in rate uncertainties of 9.4x1 O' 6 
rad/sec. This implied Q» was of the order of 10 l0 Ar rad 2 /sec where A/ is the telemetry period. A related 
analysis using the convergence properties of the Markov model resulted in Q b of the order of 2xlO' 10 Ar 
rad 2 /sec 3 . Tuning was then accomplished by: (1) choosing a constant crfrom Table 1 as stated above, (2) 
setting Q, equal to Qb during all of the runs, and (3) varying this adjustable single adjustable parameter Q lb 
about the numerical value of 10 10 A/ until the attitude errors were minimized. The accuracy of the tuning 
parameters was verified later by studying the performance of the filter over several orders of magnitude of 
Qj,b. For each dataset a few runs were also made using different a but none yielded better performance. 
Each dataset was also studied using different telemetry periods. All in all a few hundred runs were made for 
each dataset, and only a small portion of the results are shown below. 

A striking difference between the RTSF performance for RXTE and past experience with 
SAMPEX and ERBS data is the sensitivity to the numerical values of the numerical parameters, which in 
turn were somewhat dependent on the telemetry period. Thus, whereas it was sufficient for Q ib to be 
accurate to one significant figure for SAMPEX and ERBS, it turns out that the tuning parameters have to be 
accurate to three to four significant figures for RXTE. As an example, for the 11/6/97 data with 8 sec 
telemetry period, the total RTSF attitude error was 15.6 deg when Q», b = l.lxlO 10 whereas this error 
dropped to 5.1 deg when Q lb = l.OlxlO' 10 . 

TAM-Only Results 

DADMOD and RTSF (using IIC) attitude results for the 7/4/96 RXTE data with a telemetry period 
of 40 sec are presented in Fig. 1, which shows a spacecraft maneuver about the z-axis between 1500 sec and 




Figure 1. GCI-to-Body 1-2-3 Attitude Euler Angle Results for Adjusted 7/4/96 Data (solid = true, 
dashed = RTSF, crosses = DADMOD 1st root, circles = DADMOD 2nd root) 
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2500 sec. In this figure the lines represent the truth and the RTSF solutions (solid and dashed respectively) 
while the symbols represent the DADMOD correct and spurious solutions (crosses and circles respectively). 
The filter was started with IIC (large initial errors) and converged within 100 sec to a metastable spurious 
state that also shows up in the DADMOD solutions. The RTSF converges to the correct solution only about 
1000 sec later-towards the start of the maneuver. This slow convergence of the filter is a direct result of the 
inertial-pointing nature of the spacecraft, which results in the orbital motion being the sole cause for 
changes in B . (In fact, Bf is approximately constant over the maneuver period.) Note that: (1) this is the 
first independent confirmation of the DADMOD spurious solutions, and (2) such ambiguities will not arise 
if gyros provided the rate information and a TAM is used solely for attitude information. Note also a 
relatively large spread of the physical deterministic solution as a result of relatively low rate of change of 

B. 

The slow convergence severely limits any rating of the accuracy of the filter: statistics for the last 
15 points in Fig. 1 reveal r-m-s attitude errors of (0.43, 0.39, 0.17) deg about the three body axes. For more 
reliability, the filter was studied using data from the 1/4/96 dataset where an inertial span of nearly 4500 sec 
duration precedes the maneuver. The results are shown in Figs. 2 and 3. 

Fig. 2 presents sample attitude and rate results for the 1/4/96 data of 8 sec telemetry period with 
the RTSF using IIC. These results were obtained with the numerical values of Q, = Qb = 2.23x10' and 
were deemed the optimal parameters after examining the error statistics. (See Table 2 below). We see that 
the filter converges by about 4000 sec even though the initial errors ranged from about 65 deg in Angle-3 to 
about 1 13 in Angle-1. Additional residual results from the same run are presented in Fig. 3. The RTSF state 
vector evolves so as to minimize all these quantities, and we see that all are small only after 4000 sec. 

The convergence is slow here also, and it is instructive to examine the RTSF errors after 3200 sec 
separated into before, during, and after the maneuver. These are presented in the first and second columns 
of Table 2. Table 2 also compares these error statistics with the ones obtained using raw TAM data and a 
different set of tuning parameters (Q a = Qb = 1.12xl0' 9 ) separately determined to be optimal for the raw 
Hata Some clear inferences can be draw from examining Table 2. 

• The attitude errors are significant before the maneuver but noticeably decrease during the maneuver, 
which we attribute to the increased observability of changes in B . 

• In contrast to the attitude errors, the errors in the rates increase during the maneuver, which we 
attribute to (small) dynamical modeling errors of the spacecraft. 

• TAM calibration significandy improves TAM-only accuracies. 

The above information is also seen qualitatively in the bottom plot of Fig. 2 and the middle and 
bottom plots of Fig. 3. Thus, Fig. 3 clearly shows us that the rates have converged well before the attitude; 
this is in accord with past experiences with the RTSF. It should be noted, though, that there will always be 
differences in the convergence times of the RTSF attitude and rate estimates because the rates are corrected 
based on the TAM residuals. Nevertheless, it would be interesting to examine in the future if better 
accuracies result from tuning Qi and Qb to yield the same convergence times for both attitude and rates. 

The RTSF performance using IIC was further examined using the wholly inertial span of 11/6/97, 
IIC, and some of the results are shown in Fig. 4. It is clear (especially from the TAM angle plot) that the 
convergence time is of the order of 4400 sec. The error statistics from 4400 sec to the end of the data span 
are as follows: 

r-m-s attitude errors = (0.54, 0.13, 0.33) deg 

r-m-s errors in rates = (0.0049, 0.0010, 0.0024) deg/sec 
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Figure 2. RTSF Attitude and Rate Results for Adjusted 1/4/96 Data (circles = RTSF and solid = truth 

in top two plots) 
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Figure 3. Additional RTSF Results for Adjusted 1/4/96 Data Showing TAM-1 Residuals (top two 
plots) and RTSF Rate-Error Estimates (bottom plot) 
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Table 2 

RTSF ACCURACIES FOR 1/4/96 DATA SHOWING STATISTICS AFTER 3200 SEC 



Adjusted Data 

Raw Data 


R-M-S Attitude 
Errors 

R-M-S Errors in 
Rates 

R-M-S Attitude 
Errors 

R-M-S Errors in 
Rates 


(x, y, z) deg 

(x, y, z) deg/sec 

(x, y, z) deg 

(x, y, z) deg/sec 

Before 

(1.15,0.51,0.68) 

(0.0025, 0.0015, 
0.0019) 

(8.56, 2.37, 3.12) 

(0.0099, 0.0035, 
0.0075) 

During 

(0.70, 0.31,0.36) 

(0.0026, 0.0041, 
0.0037) 

(0.74, 1.09, 0.92) 

(0.0023, 0.0073, 
0.0030) 

After 

(0.38, 1.33, 0.65) 

(0.0015, 0.0038, 
0.0034) 

(0.68, 2.06, 0.57) 

(0.0017, 0.0043, 
0.0033) 
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Figure 4. RTSF Results for Adjusted 11/6/97 Data (circles = filter in top plot) 
CONCLUSIONS 

In contrast to our past experiences with SAMPEX and ERBS, which were momentum-biased 
spacecraft spinning at the orbit rate, the RTSF performance for the zero-momentum, inertial-pointing RXTE 
is characterized by extreme sensitivity to filter tuning and long convergence times (about 4000 sec). Thus, 
while past SAMPEX and ERBS results demanded accuracies of only 1 significant figure in the tuning 
parameters, it was clear that accuracies of three to four significant figures were needed for success 
application to RXTE. The performance of both DADMOD and RTSF improved significantly once the 
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telemetry period was increased from the nominal value of 2 sec to between 8 and 40 sec. Presently we 
attribute this sensitivity to a combination of telemetry period and the zero-momentum nature of RXTE. 

Careful tuning of the RTSF demonstrated per-axis attitude accuracies between 0.13 and 0.54 deg 
and rate accuracies between 0.0010 deg/sec 0.0049 deg./sec when adjusted data were used. The 
corresponding values during the maneuver of 1/4/96 were 0.31 - 0.70 deg and to 0.0026 - 0.0041 deg/sec. 
These results are similar to our past results for ERBS 8 : attitude accuracies within 0.4 deg and rate 
accuracies within 0.005 deg/sec We consider these accuracies preliminary in view of the long convergence 
times and the limited number of post-convergence results. More definitive accuracy studies are needed 
using long spans of inertial data. It would also be useful to study the relationship between the performance, 
the tuning, and the telemetry period. 

Accurate TAM calibration was performed using a recently developed algorithm, and the 
subsequent RTSF TAM-only results show a significant improvement in the attitude accuracies upon using 
calibrated TAM data. This conclusion is in general agreement with studies using a TAM-gyro 
combination 14 . 

An important theoretical development presented here is the transformation of a DADMOD 
transcendental equation into an 8th order polynomial. The first independent evidence of the accuracy of the 
DADMOD spurious solutions was obtained in the form of initial convergence of the RTSF to the 
metastable solution using the 4/7/96 data. To improve the DADMOD performance for inertial-pointing 
modes, a more sophisticated technique for smoothing time derivatives of the measured geomagnetic field 
should be applied to separate the physical root from the numerical noise. 
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